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HIGHLIGHTS 


• The effects of inner diameter and position of PA on the engine were investigated systematically. 

• PA with suitable position and inner diameter can increase pressure amplitude and acoustic power of the engine. 

• Pressure drop generated by PA and the variation of sound direct current in the loop were also analyzed. 

• The higher acoustic power is accompanied by the bigger sound direct current. 
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A thermo-acoustic Stirling engine can convert thermal energy into useful acoustic power without 
moving parts, which can be used as driving power for the refrigerators and the generators. However, the 
low efficiency of the thermo-acoustic engine limits its commercial application. In this paper, to promote 
the performance of the thermo-acoustic engine, a phase adjuster (PA) is introduced by narrowing the 
cross-sectional area on the basis of the previous studies. The influences of the position and inner 
diameter of PA on the resonant frequency, pressure amplitude, acoustic power, sound direct current (DC) 
and efficiency are studied systematically. It is found that the maximal pressure amplitude, acoustic 
power and efficiency with optimal PA could be increased by 3%, 4.5% and 10% respectively compared to 
those without PA. In addition, the effect of PA on DC flow is also discussed in detail. Results show that the 
sound DC decreases with the farther position and smaller inner diameter of PA. Therefore, the suitable 
inner diameter and position can improve the efficiency of the engine, which needs to be weighed when 
designing or optimizing PA. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

A thermo-acoustic engine based on the thermo-acoustic effect 
can realize the heat-to-sound conversion without moving parts 1 . 
In devices of this type, an acoustic wave interacts with internal 
solid structures such as regenerators and stacks to create acoustic 
power due to an imposed temperature gradient. In addition, it can 
operate with working fluids like noble gases that are benign to the 
environment. The low-grade energy such as the solar energy and 
the industrial waste heat can be used as thermal sources to produce 
the temperature gradient combined with an ambient heat 
exchanger. Therefore, developing the thermo-acoustic technology 
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is meaningful to mitigate the energy crisis and to keep the 
ecological balance. 

In a thermo-acoustic engine, the working fluid oscillates at a 
fixed frequency when the temperature gradient in the regenerator 
(stacks) reaches a critical value. In the early study of the thermo¬ 
acoustic engine, the structural size was very large, the onset tem¬ 
perature was high and the sound intensity density was usually low. 
As the new strategies and novel structures, such as pressure 
disturbance, multi-stage traveling wave thermo-acoustic engine 
and so on, were proposed, the thermo-acoustic engine has devel¬ 
oped toward the miniaturization, low onset temperature difference 
and high sound intensity. In a viewpoint of saving energy, the onset 
temperature difference should be as low as possible. Therefore, 
improving the efficiency of heat-to-sound energy conversion has 
become the subject of much research. To achieve this goal, studies 
on the thermo-acoustic effect, unsteady heat transfer in the porous 
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media, and acoustic impedance matching have been carried out by 
numerical and experimental methods [2-6 . 

In a thermo-acoustic system, the heat-to-sound conversion is 
controlled by two parameters: cat (as defined in Eq. (1)), and the 
phase angle between sound pressure and velocity [7 . 



where oj denotes angular frequency, r is the thermal relaxation 
time, r is the hydraulic radius of the regenerator, and 5 a is the 
thermal penetration thickness. For the standing wave thermo¬ 
acoustic engine, the stack plate space is about 2-4 times in thick¬ 
ness of the thermal penetration to induce a significant time delay 
between the gas motion and the gas thermal expansion (or 
contraction) to meet the Rayleigh's criterion. For the traveling wave 
thermo-acoustic engine, the hydraulic radius is much smaller than 
the thermal penetration thickness to ensure a near isothermal 
condition in the regenerator. 

The phase angle is the other parameter to determine the effi¬ 
ciency of heat-to-sound energy conversion [ 8 , 9 ]. Ceperley found 
that the efficiency of heat-to-sound energy conversion was higher 
when the phase angle between pressure and particle velocity was 
around zero 10 ]. Yazaki succeeded in realizing this thermo¬ 
acoustic setup by using the looped tube in 1998 [11 . Biwa stud¬ 
ied the thermo-acoustic phenomena in detail by combining 
experimental methods with computing methods [ 12 ]. He identified 
that the efficiency would be higher if the phase angle between 
pressure and particle velocity was closer to 0 . In consideration of 
the importance of phase angle, a number of strategies had been 
proposed to adjust the phase angle. Historically, the first study on 
the phase angle adjustment could date back to Swift [13 . In his 
research about the thermo-acoustic Stirling engine, the jet pump 
was applied to reduce the Gedeon streaming in the loop and ach¬ 
ieved a great success [14 . Nevertheless, the sudden change of the 
cross-sectional area brought by the jet pump could cause serious 
energy dissipation. Although the adjustment for phase angle was 
not mentioned explicitly, the jet pump acted as the rudiment of 
phase adjuster. Barton et al. analyzed power dissipation and time- 
average pressure through the jet pumps with different geometries 
[ 15 ]. Petculescu et al. investigated the nonlinear effect and minor 
loss of the jet pump and discussed the influence of different cone 
half-angles 16 . In recent years, Sakamoto had done much research 
on phase angle adjustment methods and proposed several feasible 
approaches for phase angle adjustment [ 17 , 18 ]. For example, one 
approach was to apply a diverging tube in the looped thermo¬ 
acoustic engine to adjust the phase angle. The temperature differ¬ 
ence of both ends in the stack reduced to 345 °C from 845 °C 19 ]. 
Another was to introduce a length of circular tube with sudden 
contraction of cross-section which is called phase adjuster (PA) in 


the thermo-acoustic cooling system. PA adjusted the phase differ¬ 
ence by narrowing the cross-sectional area, thus improving the 
cooling effect. Subsequently, they investigated the influence of in¬ 
ner diameter of PA on heat-to-sound energy conversion efficiency 
in a Loop-Tube-Type thermo-acoustic system [19 . Sahashi et al. 
studied further the working mechanism of PA in a loop-tube-type 
thermo-acoustic cooling system by experiments [20]. However, it 
was challenging to control the phase angle between pressure and 
particle velocity because of the self-excitation oscillation of the 
working fluid. To date, the working mechanism of PA was not 
completely clear. To the best of our knowledge, studies on the 
thermo-acoustic Stirling engine with PA were not reported in the 
literature. 

In this paper, inspired by the jet pump and phase adjuster, we 
introduce a PA in the thermo-acoustic Stirling engine by narrowing 
the cross-sectional area. The abrupt change of the cross-sectional 
area with PA is not larger than that with the jet pump, resulting 
in fewer energy losses than those generated by the jet jump. Firstly, 
the phase adjustment of the PA is verified. Secondly, the influences 
of the position and inner diameter of PA on the resonance fre¬ 
quency, pressure amplitude, acoustic power and the efficiency of 
the TASHE are studied systematically. The impact of PA on sound DC 
is then analyzed. Finally, the formation of vortexes is also discussed 
through the streamlines pattern to further understand and help the 
optimal design of PA structure. 

2. Numerical model 

In this section, the CFD model for the thermo-acoustic Stirling 
engine with a PA is described. The governing equations, initial and 
boundary conditions and numerical algorithm analysis of numeri¬ 
cal model are also presented. 

2.1. CFD model of the thermo-acoustic Stirling engine 

Fig. 1(a) illustrates a thermo-acoustic Stirling engine with a PA. 
The engine is comprised of a loop and a resonator, which has the 
advantages of both standing wave and traveling wave. The loop 
includes a HHX, a regenerator, a TTB and two ambient heat ex¬ 
changers. PA is placed on the top of the MAHX by narrowing the 
cross-sectional area. Inspired by the work of Lycklama [21 , an 
axisymmetric grid was applied for the modeling the system to save 
the computational time. The topological structure of the thermo¬ 
acoustic Stirling engine is performed in Fig. 1(b). The original loop 
is composed of two concentric tubes. The dimensions of the main 
parts are given in fable 1. 

Fig. 2 presents the three-dimensional schematic illustration of 
PA. L is the length of PA, and d is the diameter of PA which is smaller 
than that of the tube. Similar to the jet pump [14], the working gas 
will make a transition from a large area to a small one when flowing 



Fig. 1 . Schematic illustration of thermo-acoustic Stirling engine (a), and its two dimensional topological structure (b). 1: main ambient heat exchanger (MAHX); 2: regenerator; 3: 
hot heat exchanger (HHX); 4: thermal buffer tube (TBT); 5: secondary ambient heat exchanger (SAHX); 6: resonator; 7: feedback tube; 8: PA. 
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Table 1 

Dimensions and mesh type of main parts in the thermo-acoustic Stirling engine. 


Component 

Diameter (m) 

Length (m) 

Mesh type 

Mesh size (mm) 

Feedback tube 

0.080 

2.980 

Quad Pave 

2 

MAHX 

0.080 

0.045 

Quad Map 

1 

Regenerator 

0.080 

0.080 

Quad Map 

0.6 

HHX 

0.080 

0.058 

Quad Map 

1.2 

TBT 

0.080 

0.240 

Quad Pave 

2 

SAHX 

0.080 

0.022 

Quad Map 

1 

Resonator 

0.080-0.200 

5.000 

Quad Pave 

4 


through PA. On one hand, this transition creates a strong flow 
condition where velocity is accelerated, resulting in the variation of 
pressure and the phase between pressure and velocity. On the other 
hand, the transition will bring about an additional pressure drop 
and dissipation, called “minor losses”. To analyze the influence of 
PA on the system performance systematically, different positions 
and inner diameters of PA are investigated and compared. All po¬ 
sitions of PA are given relative to the ambient end of the regener¬ 
ator (marked by position of x in Fig. 1 ). The main dimensions of PA 
are shown in Table 2. 


2.2. Governing equations 

In a traveling wave thermo-acoustic engine, the hydraulic 
diameter in the regenerator should be much smaller than the 
thermal penetration for ensuring a good thermal contact between 
the fluid and solid. Porous zone in FLUENT satisfies this condition 
and is applied in the regenerator and heat exchangers successfully. 
The regenerator is assumed to be rigid because the pressure 
oscillation only occupies 0.03 of the mean pressure and will not 
cause the deformation of the regenerator. When the working fluid 
flows through porous zone, pressure drop can be considered by 
adding the additional momentum loss into the momentum equa¬ 
tion, which can avoid the difficulty in modeling complicated 
physical structure. The additional momentum loss can be expressed 
as [22]: 

~Si = -(y+^ c 2P f \V\v) ( 2 ) 

where S[ represents the additional momentum loss in the porous 
zone; /?,C 2 and g denote the penetration coefficient, inertial resis¬ 
tance coefficient and dynamic viscous, respectively. Some related 
coefficients applied in the porous zone are given in Table 3. 

With Eq. ( 2 ), the mass, momentum and energy conservation 
equations of the working fluid in the porous zone are given by 
Ref. [ 22 ]: 


d( £ p f ) 


1 0 


0 


at 


+ 7 Qp( er Pf v f) + ^(wx) = o 


Ox 




Table 2 

Dimensions and position of PA. 


Position — away from The ratio of inner diameter of PA to Length 

the ambient end of the the diameter of the tube (d/D) (m) 

regenerator (m) 


0.43 

0.375 

0.5 

0.625 

0.75 

0.875 

0.045 

0.52 

0.375 

0.5 

0.625 

0.75 

0.875 

0.045 

0.62 

0.375 

0.5 

0.625 

0.75 

0.875 

0.045 



—eVp + V- (cr) 



~V +\c 2Pf \v 


V 



0 

at 


( epfE { + (1 - e)p s Es) + V- (V(pf£f +?)) 
ekf + (1 — e)ks)vT + ( st ’ ~v 



where p, v, p and T are density, velocity, pressure and temperature, 
respectively; r is stress tensor; E and h are internal energy and 
enthalpy; k and e denote the thermal-conductivity coefficient and 
the porosity of the porous zone, r and x represent the radial and axis 
direction. The superscript denotes the vector, the subscript/and 
s denote the working fluid and the solid material of the porous 
zone. 

For other components, the mass, momentum and energy con¬ 
servation equations solved by FLUENT are as follows [ 22 ]: 


9(pf) 19/ x 9 , x 

-0T + ? 9?( r ^ f ) + a> v x)=o 


9 


Pfl ?) + V- (p = -Vp + V- (t) 


0 

dt 


( p f E) + V • (V (p f £ + p)) = V-[kVT + (i-7) 


( 6 ) 

(7) 

( 8 ) 


, p if 

h-- + T 

P 2 



Eqs. ( 2 )-( 9 ) will be combined with the gas state equation 
describing the total thermo-acoustic phenomena in a thermo¬ 
acoustic Stirling engine. 


2.3. Initial and boundary conditions 

In the numerical model, there are mainly two conditions to 
establish the temperature gradient: the constant heat flux and the 
constant temperature. The condition of constant heat flux means 
that the constant heat flux is input in the hot heat exchanger. Under 
this condition, the temperature gradient in the regenerator goes up 
gradually. Another condition of constant temperature means that 
the hot and ambient heat exchangers are set at constant temper¬ 
ature, and temperature gradient in the regenerator is decided by 


Table 3 

Porous materials and coefficients in the porous zone. 


Parts 

Material 

Porosity rate e 

P (m 2 ) 

C 2 ( m ^ 

Regenerator 

Stainless steel 

0.5 

1.89E-10 

4.4592E4 

HHX 

Copper 

0.675 

2.16E-7 

1.0627E3 

Ambient heat 
exchanger (AHX) 

Copper 

0.5 

2.16E-7 

1.0627E3 
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the thermodynamic process between the gas and the regenerator's 
material. In this paper, to conveniently compare the performance of 
the TASHE with/without PA and analyze the influence of PA, the 
structural parameter and the temperature of the hot and ambient 
heat exchanger are set to be unchangeable. In order to save 
computing time, the initial condition of linear temperature 
gradient in the regenerator is given originally in Table 4. Finally, the 
stable temperature distribution in the regenerator is obtained by 
solving energy equation iteratively. 

Helium gas (He) is used as the working gas. Table 4 shows some 
properties of He and initial conditions used in CFD simulation. For 
boundary conditions, all walls are assumed to have a fixed tem¬ 
perature T = 300 K except for those of the HHX, the regenerator, 
and the TBT which are set to be adiabatic. In addition, no-slip 
boundary conditions are imposed at all walls. 


2.4. Numerical algorithms analysis 


For the compressible and unsteady system in this paper, fixed 
time step with pressure-based solver and second-order implicit 
time differencing is used [23]. In our double-precision parallel 
calculation, Pressure-Implicit with Splitting of Operators (PISO) 
algorithm and second-order upwind spatial difference are selected 
to obtain the detailed information regarding the pressure and ve¬ 
locity fields. Different from the steady flow, the velocity in the 
unsteady flow is metabolic, so the Reynolds number based on the 
amplitude of velocity is used to judge the flow condition just as the 
Reynolds number in the steady flow. The Reynolds number based 
on the amplitude of velocity is defined as [24]: 


Re 


p(u)d 


( 10 ) 


where p is the density of the working gas, d is the diameter of the 
tube, /Ji is the viscosity coefficient, <u> is the velocity amplitude. 
After the system is under steady condition, the Reynolds numbers 
of the elbow and the resonator junction are 9,600 and 14,400 
respectively, indicating that the flow has been fully turbulent [24]. 
Therefore, the standard turbulence k—e model used for this model 
is valid for the fully developed turbulent flow [22]. 

To investigate the grid independence, different models with five 
kinds of meshes with 12,566, 23,900, 62,964,103,847and 128,610 
nodes are built and compared. The number of grid with 103,847 is 
chosen since the results for further refining grid has little effect on 
the results. Considering the function of different parts and their 
importance in the thermo-acoustic effect, the different mesh 
element sizes and types have been used in different parts. The mesh 
type and size is shown in fable 1. The pressure, velocity and tem¬ 
perature information can be obtained by solving governing equa¬ 
tions iteratively, so the selection of convergence criteria is very 
important. To ensure the calculation accuracy, the absolute 
convergence criteria for mass, momentum, energy, k and e equa¬ 
tion, are IE-05, IE-03, IE-06 and IE-03, respectively. 


Table 4 

The properties of gas and initial conditions used in CFD simulation. 


Physical 

Density 

Ideal gas 

properties 

Viscosity 

Constant 


Thermal conductivity 

Constant 

Initial 

Initial mean pressure 

1.5 MPa 

conditions 

Initial velocity 

0 


Temperature of the HHX 

800 K 


Temperature gradient 
in the regenerator 

300 + 6250 x (x- 1.073) 


temperature of the AHX 

300 K 


3. Pressure drop and phase modulation of phase adjuster 

3.1. Pressure drop 

To calculate the time-averaged pressure drop through PA, the 
local resistance loss coefficients I<b and K$ in the alternating current 
are assumed to be the same as that in the steady flow. This 
assumption was verified by Petculescu [16 . Although flow direc¬ 
tion varies repeatedly in one cycle, the time of forward and back¬ 
ward direction is equal. Therefore, we can divide one cycle into two 
half cycles and neglect the variations in density, which is presented 
in Fig. 3. 

According to the dynamic Bernoulli equations, during one half of 
the cycle shown in Fig. 3(a) [15 , 

a 

Pb + 2P U b-2 KspU * +p J ^ dx = Pa + 2 pU a ( 11 ) 

b 

c 

p a +^pul-^I< B pul=p c + p J -Ad x + ^ P u 2 c ( 12 ) 

a 


During the other half of the cycle as shown in Fig. 3(b), the flow 
direction changes reversely. The dynamic Bernoulli equations are 
given by: 

a 

Pc + 2 pu\ - 2 KsP u a + P j =Pc+ 2 pul (13) 

C 

b 

p a +^pul-^I< BP ul=p b + p [ALdx + r U 2 b (14) 

a 


where p a , pb, Pc Are pressure at the position of a, b, c; u a , u^, u c are 
velocity at the position of a, b, c. /C B and I <s donate the local resis¬ 
tance loss coefficients. 

Combining Eqs. (11)—(14), the pressure drop can be given as 
following: 


A p 



2 + 


( u b,max + u l 


max 


a 


(15) 


where 


1 


a 


Tu 2 

1 u max 


u 2 dt 


(16) 


o 


If the flow is sinusoidal, then a = y 2 and Eq. (15) will become 


A p 


pu 


a, max 


(/C B - I < s - 2^ + max + U c,max) 


(17) 


where ii a , m ax, Ub,max and zi c , m ax are the maximum velocity at the 
position of a, b , c. 

3.2. Phase modulation 


Sakamoto et al. [17] proposed diverging tube in a Loop-Tube- 
Type thermo-acoustic system to control the phase between pres¬ 
sure and velocity by moving the anti-node of pressure. Similar to 
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(a) (b) 

Fig. 3. Alternating current flows through PA periodically. One cycle is divided into two 
half cycles according to the flow direction, respectively shown in (a) and (b). 


the diverging tube, the introduction of PA increases the disturbance 
by narrowing the cross-sectional area, resulting in the change of 
phase between pressure and velocity. The pressure node with and 
without PA is shown in Fig. 4. The intersections of two solid lines 
and two dotted lines denote the pressure node with and without 
PA, respectively. It can be seen that the position of the pressure 
node with PA moves relative to that of the case without PA, which 
indicates the phase between pressure and velocity is modulated. 

4. Results 

4.2. Resonant frequency 

Fig. 5 shows the variation of the resonant frequency with the 
position and inner diameter of PA. The dotted line represents the 
resonant frequency 56.5 Hz without PA. It can be seen that the 
frequency change after using PA is small, less than 0.3 Hz. According 
to the formula of the resonant frequency /= a/A, the frequency is a 
function of sound velocity a and wavelength A. Sound velocity as¬ 
sociates closely with the type of gas and the system temperature. A 
is determined by the dimension of the resonator. Overall, the 
introduction of PA has little impact on resonant frequency because 
the structural change after using PA is so small that it does not 
arouse the larger change of the resonant frequency. 

4.2. Pressure amplitude 

The variation of pressure amplitude with the position and inner 
diameter of PA at the pressure anti-node is presented in Fig. 6 . The 
dotted line represents the pressure amplitude without PA. For a 
fixed position of PA expect for 0.62 m, the pressure amplitude in¬ 
creases initially and then decreases slightly. It can be seen that the 
pressure amplitude achieves its maximum value when d/D is 0.75. 
When PA is 0.62 m away from the hot end of the regenerator, 




d/D 

Fig. 5. The variation of resonant frequency with the ratio d/D at the different position 
of PA. 


pressure amplitude increases monotonously with inner diameter of 
PA, but the increase becomes flat gradually. For a fixed inner 
diameter of PA, pressure amplitude decreases with an increase in 
the position of PA except PA with d/D = 0.875. Nevertheless, this 
effect weakens as inner diameter of PA increases. It is observed that 
the difference of pressure amplitude at three positions of PA is 
small when the d/D is 0.875. Overall, the maximal pressure 
amplitude increases by 3% compared to the case without PA. 
However, pressure amplitude with the d/D = 0.375 is smaller than 
that without PA regardless of the position of PA. 

4.3. Acoustic power 

Acoustic power characterizes the capacity of heat-to-sound 
energy conversion at the constant temperature difference be¬ 
tween the two ends of the regenerator. The higher the acoustic 
power is, the stronger the thermo-acoustic effect is. Fig. 7 shows 
the variation of acoustic power with the position and inner diam¬ 
eter of PA at the cold (a) and hot (b) ends of the regenerator. The 
variation trend of acoustic power with the position and inner 
diameter of PA is similar to that of pressure amplitude. Specifically, 
the acoustic power is higher for the larger inner diameter and the 
nearer distance from the hot end of the regenerator. When the d/D 



Fig. 4. The pressure node with and without PA. 


Fig. 6. Pressure amplitude as a function of the ratio d/D at different PA positions. 
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(a) (b) 

Fig. 7. The variation of acoustic power at the cold end (a) and the hot end (b) of the regenerator with the ratio d/D. The dotted line represents acoustic power at two ends of the 
regenerator of the engine without PA. 


is 0.75, and PA position is 0.43 m from the hot end of the regen¬ 
erator, the acoustic power achieves its maximum value which in¬ 
creases by 4.5% relative to the case without PA. 


4.4. Sound direct current (DC) 


The loop in the thermo-acoustic Stirling engine generates “short 
circuit” for acoustic DC (Gedeon Current) from the hot heat 
exchanger to the secondary cold heat exchanger. Following Swift 
[13], the sound DC is given by 



Re 


PiL / 1 



2 + Pm ^20 


( 18 ) 


It can be clearly seen that the pressure drop is negative, 
implying that the U 20 brought by pressure drop can offset a part of 
acoustic DC. The absolute value of pressure drop increases with a 
decrease in inner diameter of PA, and increases with position of PA. 
According to Eq. (17), pressure drop is determined by the square of 
particle velocity and increases with particle velocity. On one hand, 
the particle velocity increases as inner diameter of PA decreases. On 
the other hand, particle velocity increases with the distance from 
the hot end of the regenerator. Obviously, the higher negative 
pressure drop is more helpful to suppress the acoustic DC. There¬ 
fore, the higher pressure drop brings about the smaller sound DC 
comparing Figs. 8 and 9. 


where M 2 denotes the sound DC; p\ and p m are the oscillating 
density and the average density; U 1 and U 2 0 are the oscillating 
velocity and the second order time-average volumetric velocity; Re 
means the real part; the tilde denotes complex conjugation. The 
first term on the right hand of Eq. (18) is nonzero as long as the 
system oscillates. If none measures are taken to produce U 2 o to 
offset sound DC, heat will flow directly from the hot exchanger to 
the secondary cold heat exchanger, causing unnecessary heat 
losses. 

Fig. 8 illustrates the variation of sound DC with the position and 
inner diameter of PA. The dotted line represents sound DC without 
PA and its value is 0.039 kg s -1 . The sound DC decreases with the 
increase in position of PA at a fixed inner diameter of PA. However, 
the influence of position of PA becomes weak as the inner diameter 
of PA increases. At a fixed position of PA, the sound DC decreases 
with inner diameter of PA. When the PA position is 0.63 m away 
from the hot end of the regenerator, the sound DC with d/D = 0.375 
decreases by 30% relative to the case without PA. Nevertheless, the 
sound DC is 0.040 kg s -1 higher than that without PA when the d/D 
is 0.75 or 0.875. Comparing Figs. 7 and 8, it can be found that the 
sound DC increases while acoustic power increases. 

When the alternating current flows through PA, additional 
pressure drop is generated from a large cross-sectional area to a 
small one (or the other way around). We use Eq. (17) to calculate 
the pressure drop, and the variation of pressure drop as inner 
diameter and position of PA is shown in Fig. 9. 


4.5. Efficiency 

The variation of the efficiency of the engine with inner diameter 
and position of PA is shown in Fig. 10. It can be seen that the effi¬ 
ciency firstly increases and then decreases with the increase in 
inner diameter of PA, and decreases with the increase in position of 
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Fig. 9. The variation of pressure drop with the ratio d/D. 



Fig. 11. Pressure wave when the oscillation stabilizes, d\, d 2 , 63 , 0 4 are chosen repre 
senting different phases in one cycle. 


PA. The efficiency of the engine with PA increases compared to the 
case without PA when the d/D > 0.5. Especially, the maximal effi¬ 
ciency increases by 10% with the d/D = 0.75, PA position of 0.43 m 
away from the ambient end of the regenerator. 

4.6. Oscillatory flow field 

Fig. 11 shows the diagram of the pressure wave while the engine 
levels off. Four points in one cycle are chosen to illustrate the effects 
at phases of d\, 62 , O 3 and 64 . The streamlines at the position of PA, 
the elbow and the resonator junction are illustrated in Fig. 12. 
When the alternating current flows through the cross-sectional 
area with a sudden change and turns around, vortexes with 
different intensity and sizes form. At the elbow, the periodic sep¬ 
aration and converge of the alternating current because of turning 
around control the formation of vortexes. The sudden variation of 
velocity due to abrupt change of the cross-sectional area intensifies 
the formation of vortexes. The number and intensity of vortexes 
decrease when the working fluid flows anticlockwise. The intensity 
and size of vortexes around PA vary with time because the alter¬ 
nating flow experiences the sudden expansion and shrink of the 
cross-sectional area back and forth. The fluid separates and con¬ 
verges periodically at the resonator junction, thus generating vor¬ 
texes. It should be noted that PA aggravates the formation of 
vortexes by narrowing the cross-sectional area. The formation and 



d/D 

Fig. 10. The variation of efficiency as the function of the ratio d/D. 


shedding of vortexes could affect the distribution of velocity and 
temperature, causing unnecessary energy losses [25,26 . Therefore, 
it is helpful to optimize the structure of PA by analyzing the cause 
and position of vertexes formation. 

Based on the previous analysis, the acoustic power decreases 
with smaller inner diameter of PA while increases with larger inner 
diameter of PA. The sound DC is smaller with smaller inner diam¬ 
eter of PA while is larger with larger inner diameter of PA. The 
reason that causes the variation of acoustic power and sound DC 
will be analyzed in the following section. 


5. Discussions 


5.1. Influence analysis of PA 


The acoustic impedance at the ambient end of the regenerator 
determines the acoustic power flow into the regenerator. And the 
acoustic power flow into the regenerator can be approximately 
expressed as 6] 



1 

2 


Re 


P1U1 



2R e[Z] 


(19) 


where Z = P 1 /D 1 is the acoustic impedance at the ambient end of 
the regenerator, p\ and Lb are pressure amplitude and volume flow 
rate amplitude. For given pressure amplitude at the ambient end of 
the regenerator, high acoustic impedance will reduce acoustic po¬ 
wer flow into the regenerator. PA can accelerate velocity by nar¬ 
rowing the cross-section area, thus affecting the volume flow rate. 
The variation of volume flow rate at the ambient end of the 
regenerator as inner diameter and position of PA is shown in Fig. 13. 

It can be seen that volume flow rate firstly increases and then 
decreases with the increase in inner diameter of PA. When the 
inner diameter of PA is fixed, volume flow rate decreases with the 
increase in position of PA. According to Eq. (19), the increase of 
volume flow rate decreases acoustic impedance at the ambient end 
of the regenerator, and then decreases the acoustic impedance. As 
such acoustic power flow into the regenerator is increased. When 
the inner diameter is smaller (d/D = 0.375), the volume flow rate 
decreases compared to the case without PA, thus acoustic power 
flow into the regenerator is reduced. PA with the larger inner 
diameter (d/D > 0.6) increases the acoustic power flow into the 
regenerator. 















1146 


Y.-W. Liu, P. Yang / Applied Thermal Engineering 73 (2014) 1139—1148 




Fig. 12. Streamline at different moments in one cycle for three different positions: the elbow (left), PA (middle), and the resonator junction (right). 


However, when the alternating current flows through PA, it will 
make a transition from a large area to a small one, which can 
produce additional pressure drop and power dissipation called 
“minor losses”. The minor losses can be evaluated by a resistance, 
which is given by 14] 



Fig. 13. The variation of volume flow rate at the ambient end of the regenerator with 
the ratio d/D. 


< 20 > 

where R P a denotes the acoustic resistance of PA, as represents the 
cross-sectional area of PA, p m and u is the mean density and velocity 
of working gas at the position of PA, JC B and I < s are local resistance 
loss coefficients. Fig. 14 shows the variation of the resistance with 
the position and inner diameter of PA. 

It can be seen that the resistance decreases with the increase in 
inner diameter of PA, and decreases with the position of PA. 
However, the impact of the position of PA weakens as inner 
diameter of PA increases. The power dissipation increases with R PA . 
Therefore, the change of the resistance replacing PA reflects the 
variation of power dissipated by PA. It is not surprising that the 
smaller resistance means less power dissipated by PA, which has 
less effect on the engine. 

Overall, the acoustic power is determined by the two combined 
effects above. On one hand, PA changes the acoustic impedance at 
the ambient end of the regenerator by narrowing the cross- 
sectional area. On the other hand, PA also dissipates power which 
can be measured by PA's resistance. When the inner diameter of PA 
is smaller (d/D = 0.375), the decrease of volume flow rate results in 
the decrease of acoustic power flow into the regenerator, and at the 
same time power dissipated by PA is maximal. The two effects both 
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Fig. 14. The variation of resistance of PA with the ratio d/D. 


worsen the performance and result in the obvious decrease of 
acoustic power. When inner diameter of PA is larger, the volume 
flow rate increases and power dissipated by PA decreases, which 
implies the more acoustic power flow in the regenerator gets the 
run upon the power dissipated by PA. Therefore, the performance is 
promoted. In addition, the volume flow rate increases and Rpa de¬ 
creases with the decrease of PA's position. The closer position away 
from the hot end of the regenerator is more helpful to improve the 
performance. In summary, the larger inner diameter and closer 
position of PA is favorable for the system. 


acoustic impedance at the ambient end of the regenerator is 
changed, which determines the acoustic power flow into the 
regenerator. On the other hand, PA also dissipates part power 
because of the sudden change of the cross-sectional area. The two 
combined factors both determine the effect of PA on the perfor¬ 
mance. In this study, the impact of these factors was investigated 
systematically by changing inner diameter and position of PA. The 
conclusions can be drawn as follows. 

(1) The variation of resonant frequency with PA is small despite 
how the position or inner diameter of PA changes. The reason 
is that the resonant frequency is mainly determined by the 
resonator, and the structural change after using PA is so small 
that it does not arouse the larger change of the resonant 
frequency. 

(2) Pressure amplitude and acoustic power are higher with the 
larger inner diameter of PA and the nearer distance from the 
hot end of the regenerator. Larger inner diameter and closer 
position of PA decrease the acoustic impedance at the cold of 
the regenerator, which increases the acoustic power flow 
into the regenerator. Concurrently, PA with larger diameter 
and closer position dissipates less power. Thus, the perfor¬ 
mance is promoted. The maximum values of pressure 
amplitude, acoustic power and efficiency increase by 3%, 4.5% 
and 10%, respectively. 

(3) Sound DC in the loop decreases with inner diameter of PA, 
but pressure drop generated by PA increases with inner 
diameter of PA. The bigger pressure drop can suppress the 
sound DC more effectively. When the inner diameter of PA is 
30 mm, acoustic DC decreases by 28%. 


5.2. Analysis of sound DC 

Sound DC in the loop inevitably causes unnecessary heat losses, 
lowering the efficiency of heat-to-sound energy conversion. To 
suppress sound DC is an important way to improve the perfor¬ 
mance of thermo-acoustic Stirling engine. 

At the cold end of the regenerator, we can rewrite the first term 
on the right hand of Eq. (18) following Gedeon [27] 
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P\U\ 




So sound DC can be rewritten 


( 21 ) 


^2 — PmWc/pm + Pm^20- (22) 

According to Eq. (22), the variation of sound DC is closely related 
to acoustic power W c and the second order volume flow rate U 20 . 
On one hand, sound DC increases with acoustic power. On the other 
hand, the alternating current flowing through PA can generate 
pressure drop, producing U 20 which can counteract a part of sound 
DC. When inner diameter of PA is smaller, acoustic power is mini¬ 
mal and pressure drop is maximal, leading to remarkably small 
sound DC. Sound DC with the ration d/D = 0.5 is determined by 
relative change of W c and U 2 0 . If the increase of sound DC caused by 
W c is smaller than the reduction of sound DC offset by U 20 , sound 
DC would decrease. Sound DC with PA of the ratio d/D = 0.625, 
0.75and 0.875 increases mainly because of the increase in W c . 
Under these conditions, pressure drop produced by PA is so small 
that U 20 cannot offset the increase in sound DC. 


6. Conclusions 


PA creates a strong flow condition in which particle velocity is 
accelerated by narrowing the cross-sectional area. On one hand, the 


Nomenclature 


Symbols 

C 2 inertial resistance coefficient 

d diameter of PA, mm 

D diameter of tube, mm 

E internal energy, W 

h enthalpy, W 

R resistance, kg m~ 4 s _1 

k thermal conductivity coefficient, W m” 1 k 1 

L length, m 

a s cross-sectional area of PA, m 2 

I< local resistance loss coefficient 

p pressure, Pa 

r hydraulic radius, mm 

Si source, kg s~ 2 

T temperature, K 

u particle velocity, m s -1 

U volume flow rate, m 3 s -1 

x axis position, m 

Ap pressure drop, Pa 

a sound velocity, m s _1 

W acoustic power, W 
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Greek symbols 

16 penetration coefficient 

ju coefficient of kinetic viscosity, kg m -1 s -1 

r stress, Pa 

co angular frequency, s -1 

v velocity, m s -1 

p density, kg m -3 

(p phase angle, ° 

Superscripts 

vector 

complex conjugation 
Subscripts 

c cold end of the regenerator 

f working fluid 

x axial direction 

res resonator 

B expansion 

m mean 

h hot end of the regenerator 

s solid 

r radial direction 

reg regenerator 

S contraction 

Other 

|| fluctuation amplitude 
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